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Abstract. In the last decade, numerous Virtual Observatory organizations were es- 
tablished. One of these is the German Astrophysical Virtual Observatory (GAVO) that 
e.g. provides access to spectral energy distributions via the service TheoSSA. In a pilot 
phase, these are based on the Tubingen NLTE Model-Atmosphere Package (TMAP) and 
suitable for hot, compact stars. We demonstrate the power of TheoSSA in an application 
to the sdOB primary of AA Doradus by comparison with a "classical" spectral analysis. 



1. Introduction 

With the progress in observational and computer technology, the quantity of astronom- 
ical data was and is increasing exponentially. To be able to handle the actual and future 
amount of data, Virtual Observatories (VOs) were established. Their aims are to make 
the data globally accessible, to prevent losing data, and to develop tools, services, and 
standards to enable research using the internet. The idea of the VO is rather old, start- 
ing with the first databases and tools such as e.g. SIMBAE0(1979), ALADIN 1 (1999), 
or IRAF 1 (1984). Since 2000 more and more national VOs were founded to make as- 
tronomical data accessible. Many of them joined the International Virtual Observatory 
Alliancd3 (IVOA). Observational and theoretical data sets were prevented from getting 
lost by making them public. Standards were defined to get a uniform shape of all data. 
Services and tools for data-handling were developed, and contact points were created 
to aid publishing data and using tools and services. 

One of these national organizations is the German Astrophysical Virtual Obser- 
vatory 

(GAVG0), 

the German effort to contribute to this idea. Within the framework 
of a GAVO pr oject, the well- established Tubingen NLTE Model-Atmosphere Package 
(TMAPE Werner et al. 2003) was made accessible. This model-atmosphere code and 



the different ways to access its products are described in the following sections. 



'http://simbad.u-strasbg.fr/simbad/, http://aladin.u-strasbg.fr/aladin.gml http://iraf.noao.edu/ 

2 http://www.ivoa.net/ 

3 http://g-vo.org/ 

4 http://as tro.uni-tuebingen.de/~TMAP 
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2. The Basis - The Tubingen NLTE Model-Atmosphere Package 

TMAP is used to calculate stellar model atmospheres in hydrostatic and radiative equi- 
librium. It assumes plane-parallel geometry due to the relatively thin atmosphere of a 
hot, compact star compared to its radius. Up to about 1500 atomic levels can be treated 
in NLTE, about 4000 lines for classical model atoms and more than 200 millions of 
lines for iron-group elements (Ca-Ni) can be consid ered. Currently all elements from 
hydrogen to nick el can be consider ed dRauchl 12003). For the iron-group elements the 



program IrOnlc (IRauch & D eetien 2003) is used to calculate sampled cross-sections 
and to set up atomic data files. This program is necessary to handle the hundreds of 
millions of lines of iron-group elements and to reduce them with a statistical treatment 
to a number suitable for the model atmosphere code. 

TMAP is used for hot, compact objects with T eS > 20000 K and 4 < g < 9. The 
GAVO service TheoSSA is based on this model-atmosphere code at the moment. 



3. Implementation - The VO Service TheoSSA 

TheoSSA^ was originally planned to provide spectral energy distributions (SEDs) as 
input for photoionization models for planetary nebulae (PNe). For these models, black- 
bodies were used very often and, thus, such PN analyses yielded wrong results. The 
blackbody approximation differs vastly from our model fluxes (Fig.[TJ). In order to sup- 
ply reliable ready-to-use SEDs this service gives access in three ways. 

1) download already calculated SEDs 

2) calculate individual SEDs with the simulation software 

3) download ready-to-use atomic data or tailor individual model atoms with TMAT^\ 

The TheoSSA database is controlled via a web interface where the fundamental 
parameters (T e ff, logg, abundances) are entered. Then all available SEDs within a cho- 
sen parameter range are displayed and can be downloaded. The database is growing in 
time, because SEDs calculated with TMAW are automatically ingested. In case that the 
calculation procedure is significantly changed, these SEDs are re-calculated with the 
latest TMAW version. 



4. More Detailed/Individual SEDs - TMAW 

TMAW is the simulation software of the model-atmosphere code. Like TheoSSA, it 
is controlled via web interface where the requested parameters and the email address 
are entered. The sending process starts the calculation of a SED with the requested 
parameters on the PC cluster of the Institute for Astronomy and Astrophysics Tubingen 
or, if the number of requests is > 150, on compute resources of AstroGrid-T^. 



5 http://dc.g- vo.de/theossa 
6 http://astro.uni-tuebingen.de/~TMAW 
7 http://as tro.uni-tuebingen.de/~TMAD 
8 http://www.gac-grid.net/ 
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Figure 1. Comparison of the TMAP (thin line, blue in the online version) and 
TMAW (thick, red) fluxes (T e g= 120000 K, log g = 7.5). The strongest deviations 
can be seen in this region, but still they are very small compared to a blackbody flux 
(black) with the same T (insert). 



In the standard procedure, the steadily updated, ready-to-use model atoms from 
TMAD are taken as input and, after the LTE start model is created, a continuum model 
is calculated (Fig.©. Then lines are inserted i nitially withou t temperature correction, 
then with Unsold-Lucy temperature correction (Uns6ld~l ll968l) . After about one or two 
days and at least 300 iterations, line profiles are calculated and the result is sent to the 
email address of the user. 

At the moment, this service can consider the elements hydrogen, helium, carbon, 
nitrogen, and oxygen, but in the future the next most abundant elements, neon and 
magnesium (as well as the iron-group elements, Ca - Ni), will be included. 



5. Creating The Model Atoms - TMAD 

The Tubingen Model Atom Database {TMAD) was originally created for TMAP and 
therefore its format is compatible with it. It is updated regularly. For example the He i 
model atom was enlarged in the latest update due to recently published atomic data 
retrieved from NIST. Now levels up to a principal quantum number of n = 10 instead 
of n = 5 can be considered as NLTE levels (Fig. [3]). 
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Figure 2. Standard procedure of a TMAW calculation. The numbers are taken 
from the example given in the text. 
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Figure 3. Grotrian diagram of HEi (as provided by TMAD), singlet (left) and 
triplet (right) systems. 



Presently, the elements H, He, C, N, O, F, Ne, Na, Mg, Si, S, Ar, and Ca are 
included. Main data sources for TMAD are the standard atomic databases like 
or the Opacity Project 

TMAD supplies two kinds of data. Level energies and radiative as well as col- 
lisional transitions can be found for model atmosphere and for SED calculations (the 
latter including fine-structure splitting). For model atmosphere calculations ready-to- 
use or individual data can be downloaded. This data is accessible to the community and 
can be used for all desired purposes, e.g. as input for another model-atmosphere code. 



9 http://www.nist.gov/pml/data/asd.cfm 
10 http://cdsweb.u-strasbg.fr/topbase/op.html 
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Figure 4. Comparison of the number of POS (left) and LIN (right) lines for Fe 
(top) and Ni (bottom) before and after the extension of Kurucz' line lists. 



6. Handling Iron-Group Elements - IrOnlc 



Another GAVO project is the parallelization of the IrOnlc code dRauch & Deetjeiill2003 ). 



Presently, it may take one to three days to calculate cross-sections and a model atom 
with IrOnlc. To accelerate, we write a new version of this code using MPI as well as 
GPU (CUDA) parallelization. This service will be controlled via a web interface in a 
similar way to TMAW and will be accessible to the public. 

Input for the IrOnlc code are Kurucz' line list^] These contain LIN and POS 
lines. The latter are only lines with measured (good) wavelengths, LIN lines have 
theoretically calculated entries additionally and therefore a bout 300 time s more lines 
than POS lists. New data became available recently (Figj4j Kumc zj|2009l). Additional 



data s ources are the OPACITY and IRON projects (Seato net al.lll994l : iHummer et al. 



1993) 



IrOnlc divides the energy range up to the ionization energy of an ion into several 
bands (Fig. [5). All levels within one band contribute to this so-called superlevel via 
defining its energy and statistical weight. This superlevels are treated in NLTE, the 
levels within one band are in LTE relation. Therefore the number of NLTE levels is 
reduced to a number manageable for the model atmosphere code. For all transitions 
within one band or between different bands, the cross-sections are sampled (Fig. [6]). 



"http://kurucz.harvard.edu/atoms.html 
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Figure 5. Energy range of Fevn divided into bands to create superlevels. The 
ionization energy is 1.008 • 10 6 cirT 2 . 



7. A First Benchmark Test - The Case of AA Doradus 

The initial aim of TheoSSA and TMAW was to provide easy access to synthetic stellar 
fluxes (Sect.O for the PNe community, which was interested in realistic stellar ionizing 
fluxes for their PN models. It turned out rather quickly, that other groups were interested 
in X-ray, optical, and infrared fluxes as well. TheoSSA / TMAW were then extended 
for this purpose. While the accuracy in ionizing fluxes is better than 10 % compared 
to TMAP (Sect.O calculations with the most detailed model atoms, the precision in 
the optical is limited due to the smaller, standard model atoms that are used in the 
TMAW model-atmosphere calculations. The aim is here to perform preliminary spectral 
analyses with TMAW SEDs and to be better than 20% in T e s, logg, and abundance 
determinations. 

Publishing data of all kind is a very important task of the VO. For all these data, 
quality control should be a pre-requisite and is highly desirable. Therefore it is neces- 
sary to establish benchmark tests for VO data, services, and tools. This is essential for 
theoretical data and simulations as well. 



7.1. A A Dor 

The eclipsing, post common-envelope binary AADor has an sdOB star primary star 
whose optical spectrum shows the broad Balmer lines and the He n A 46 86A line, typ- 
ical for this class. A recent analysis of the primary ( Klepp & Rauchl [201 lh yielded 
T e f\- = 42000 K and log g = 5.46 within small error limits. For this analysis, an extended 
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Figure 6. Sampled RBB cross-sections for Fe vn (example transition between our 
bands 2 and 5, Fig . |5j . 



grid of NLTE model atmospheres considering opacities of the elements H, He, C, N, O, 
Mg, Si, P, S, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, and Ni was calculated. T e s and logg were 
then determined by a^ 2 -fit method based on synthetic spectra of these models and high- 
resolution and high-S/N optical UVES observations (cf. Rauch these proceedings). 

AADor is, thus, an ideal object to prove the TMAW service for its adequacy 
for sp ectral analysis. For this test, we recalculated the fine grid of iKlepp & Rauch 
(1201 lh via TMAW (Fig.©. We used their abundances, T e ff = 39500 - 43500 K with 
A7/ eff = 250 K, and log g = 5.30 - 5.60, with Alogg = 0.01. In total, 527 HHeCNO- 
composed models were calculated in one week on compute resources of AstroGrid- 
D. Analogously to Klepp & Rauch] (1201 lh . we calculated synthetic spectra from these 
models and performed &x 2 fit to the UVES0 observations (Progld 66.D-180). Aston- 
ishingly, the result is very good (Fi g . U}- We arrive at T^g = 41 150 K and log g = 5.43 
while the x 2 fit of IKlepp & Rauchl (1201 lh results in T s g = 40600 K and log g = 5.46. 
This are deviations of about 1 % in r e ff and 1 % in log g, which is better than our ex- 
pectation. The wider and irregular contour lines of the TMAW-based x 1 fit (Fig. [7]) 
are due to the fact that onl y about 300 iterations were performed in contrast to about 
8000/model in the grid of IKlepp & Rauchl d201 lh and, thus, the TMAW models are not 
homogeneously converged. 



Ultraviolet and Visual Echelle Spectrograph at ESO's VLT 
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7.2. Conclusion and further improvements 

The example of AA Dor shows clearly, that TMAW models are well suited for the spec- 
tral analysis of optical spectra. The calculation times for grids via TMAW are small 
because of the AstroGrid-D computational power in the background. However, with 
some improvements, we enhanced the reliability of TMAW models. 

To ensure a reasonable calculation time, TMAW uses adjusted, relatively small 
model atoms (e.g. 207 NLTE leve ls and 498 lines, Sect.H]) and the Unsold-Lucy tem- 



perature correction (lUnsold II1968I) . This temperature correction is numerically stable, 
but the structure differs from one calculated with standard corrections (Fig. [8). Al- 
though both structures mainly agree in the line-forming regions, the deviations are not 
negligible. At log/n>0, the temperature structures differ more strongly (Fig. [8) due to 
the iron-group opacities that are considered in the TMAP models. Iron-group opacities 
will be considered by TMAW in the near future. 

To compensate this deficiency, we included additional calculation steps, first using 
the standard temperature correction applied by TMAP and then increasing the number 
of atomic levels in a subsequent line-formation step by using far more detailed model 
atoms (402 NLTE levels and 1652 lines in this example). The result is shown in Fig.[9] 
With these improvements, the accuracy is improved and the calculation time rises only 
slightly. 
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Figure 8. Comparison of the TMAP (thin, blue) and TMAW (thick, red) tempera- 
ture structure (T e jf = 42000 K, log g = 5.46). The black line shows the structure after 
the post-calculation steps. The formation depths of the line cores of hydrogen and 
helium are marked. 
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Figure 9. Comparison of theoretical lift (left) and H<5 (right) line profiles cal- 
culated from TMAP (thin, blue) and TMAP (thick, red) models (T eff = 42000 K, 
log g- 5.46). HjS (left) and H<5 (right) are displayed because they are formed at dif- 
ferent depths of the atmosphere (Fig.[8j. The TMAP model shows additional metal 
lines due additional species included. A rotational velocity of 30 km/sec is applied. 
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